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Abstract. We have formulated the momentum equation 
for sectorial line driven winds from rotating stars includ- 
ing: (a) the oblateness of the star, (b) gravity darkening 
(von Zeipel effect), (c) conservation of angular momen- 
tum, (d) line driving specified by the force multiplier pa- 
rameters (k, a, S), (e) finite disk correction factors for an 
oblate star with gravity darkening for both the continuum 
and the line driving. The equations are solved numerically. 
We calculated the distribution of the mass flux and the 
wind velocity from the pole to the equator for the winds 
of B[e]-supergiants. Rotation decreases the terminal ve- 
locity in the equatorial region but hardly affects the wind 
velocity from the poles; it enhances the mass flux from the 
poles while the mass flux from the equator remains nearly 
the same. These effects increase with increasing rotation 
rates. 

We also calculated models with a bi-stability jump 
around 25 000 K, using force multipliers recently calcu- 
lated with a Monte Carlo technique. In this case the mass 
flux increases drastically from the pole to the equator and 
the terminal velocity decreases drastically from pole to 
equator. This produces a density contrast in the wind 
p(equator)/p(pole) of about a factor 10 independent of 
the rotation rate of the star. We suggest that the observed 
density contrast of a factor ~ 10 2 of the disks of B[e] stars 
may be reached by taking into account the wind compres- 
sion due to the trajectories of the gas above the critical 
point towards the equatorial plane. 

Key words: stars: early type - stars: emission line, Be - 
stars: B[e] - stars: mass loss - stars: rotation - supergiants 
- stars: winds, outflows 



1. Introduction 

In this paper we study the effects of rotation on the radia- 
tion driven winds of early-type supergiants. We will focus 



on the explanation for the occurrence of disks around fast 
rotating B[e] supergiants. 

B [c] sup ergiants, also designated sgB[e] stars (Lamers 
et al. 1998 ), are B type supergiants that exhibit forbidden 
emission lines in their optical spectra. The observations 
of hybrid spectra, i.e. spectra with broad UV P Cygni 
features and narrow emission lines and dust emission of 
B[e] supergiants in the Magellanic Clouds led to the pro- 
posal of a disk wind model by Zickgraf et al.( |l985| ). This 
model postulates a dense disk of outflowing material in a 
fast line driven wind to explain the observed character- 
istics of the spectra of these stars (see Zickgraf 1992 ). 
The stellar wind at the equator is about ten times slower 
than that at the pole. Also, the wind at the equator is 
about a hundred times denser than at the pole. Addi- 
tional evidence for a two-component outflow has been ob- 
tained from polarimetric measurements, e.g. by Zickgraf 



& Schulte-Ladbeck (1989) 



The precise mechanism behind the formation of these 
disks is still a mystery (Cassinelli 1998). It is however clear 
that the origin of an axisymmetric wind structure such as 
a disk may well be connected to the fast rotation of a star. 
Two theories using rotation in a different way have been 
considered for the formation of these disk winds: 

- (1) The Wind Co mpres sed Disk (WCD) model of 
Bjorkman & Cassinelli ( 1993| ), that invokes the kinemat- 
ics of the winds from rotating stars. The streamlines of the 
gas in the wind from both hemispheres of a rapidly rotat- 
ing star cross in the equatorial plane. The concentration of 
the gas and the shock in the equatorial plane produce an 
outflowing equatorial disk, with a thickness on the order 
of a few degrees. Owocki & Cranmer (1994) have argued 
that the motion to the equatorial plane may be counter- 
acted by the radiation force perpendicular from the plane 
due to lines. However, this "wind-inhibition effect" may 
not be effective in winds with a strong density gradient in 
the equatorial direction. We return to this in Sect. 7. 

- (2) The rotationally induced bi-stability model (RIB) 



of Lamers & Pauldrach (1991) invokes an increase in the 
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mass flux from the equator and a decrease in the equatorial 
wind velocity compared to the poles from the bi-stability 
jump. This jump in mass flux and in wind velocity is due 
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to the temperature difference between the pole and the 
equator of a fast rotating B[e] supergiant with gravity 
darkening. The jump will occur for stars with effective 
temperatures between 20 000 and 30 000 K. 

In reality, both effects, i.e. the wind compression and 
the rotation induced bi-stability, may be operating to- 
gether and amplifying one another in the line driven winds 
of rapidly rotating early-B stars. (For a detailed explana- 
tion of both models, see Lamers & Cassinelli (1999) (here- 
after ISW), Chapter 11.) 

The theory of radiation driven winds for non-rotating 
stars was developed by Castor et al. ( 1975| ) (hereafter 
CAK) and predicts the mass loss rate and Voo for spheri- 
cal winds. The influence of rotation on line driven winds 
was investigated by Friend & Abbott ( 198C| ). They found 
an increase in the mass-loss rate at the equator and a 
corresponding decrease in Voo, but these authors did not 
consider the effects of gravity darkening and the oblate 
shape of the star. Cranmer & Owocki ( |1994[ ) described 
the finite disk correction factor in case of an oblate star. 

In this paper we will modify the line driven wind the- 
ory, including these effects of oblateness and gravity dark- 
ening as well as the rotational terms in the equation of 
motion. This will then be applied to rapidly rotating B[e] 
supergiants to investigate whether these effects can ex- 
plain the occurrence of disks. We will also apply recent 
bi-stability calculations to study the effect of the RIB- 
modcl to explain the occurrence of outflowing disks around 
rapidly rotating B[e] supergiants. 

In Sect. 2 we describe the theoretical background of 
the radiation driven wind theory and the bi-stability ef- 
fect. In Sect. 3 we derive the equations for the winds from 
oblate rotating stars with a temperature gradient between 
poles and equator due to gravity darkening. In Sect. 4 we 
describe the method for solving the equations and for cal- 
culating the winds from rotating stars. These calculations 
will be applied to B[e] supergiants in Sect. ||. In Sect. || we 
will investigate the effect of bi-stability on rotating stars 
and Sect. [7| concludes with a summary of this work and a 
discussion on the formation of disks of B[e] supergiants. 



2. Theoretical context 

The computation of the dynamics of line driven stellar 
winds is a complicated problem in radiation hydrodynam- 
ics. It involves the simultaneous solution of the equations 
of motion, the rate equations and the radiative transfer 
equations in order to calculate the radiative acceleration. 

For hot luminous stars the line forces are the most 
important driving forces in the wind. CAK have shown 
that the line acceleration can be parameterized in terms 
of the optical depth parameter t ~ p(dr/dv) as g\i nc ~ kt a , 
where k and a are parameters that depend on composition 
and temperature of the wind. In this expression k is a 
measure of the number of lines and a is a measure of the 



distribution of the line strengths with a = or 1 for a 
pure mix of optically thin or thick lines respectively. 

In this paper we adopt the CAK formalism and sim- 
plify the equations of motion by assuming a stationary, 
radial flow of a viscousless fluid. The possible influence 
of magnetic fields will be ignored. These simplifications 
are subject to the following restrictions (see e.g. Abbott 
1980): 

— The CAK theory treats the line acceleration in a sim- 
plified way in terms of the force multipliers, described 
above. Detailed calculations of the line driving, includ- 
ing millions of lines, show that the total line accelera- 
tion can be well described by this simple representation 
(Pauldrach et al. 1986). 

— Both observations and calculations show that line 
driven stellar winds are not stationary. Even a wind 
in a stationary solution will develop shocks (see e.g. 
Lucy 1982 ). However, the time-averaged structure fol- 
lows the stationary state quite well (see Owocki et 
al. 1988, and Feldmeier 1999). Therefore we will re- 



strict ourselves to stationary models. 

— A sectorial model will be adopted. This means that 
for every latitude a one dimensional problem will be 
solved. Wind compression as in the WCD model will 
thus be neglected, but the main effect of wind compres- 
sion is expected to redistribute the mass loss and not 
to change the total mass that is lost from the rotating 
star. 

— The absence of viscosity is a good approximation at 
the high density of line driven winds (see CAK). 

The representation of the line force by a simple power 
law may seem to be a gross simplification of the underly- 
ing physics of myriad line absorption processes, but it can 
be shown to hold for a homogeneousl y dist ributed mix o f 
optically thick and thin lines (Abbott 1982; , Gayley |l995| ) . 
Calcula tions of CAK and more recently in NLTE by Vink 
et al. ( 1999 ) of realistic model atmospheres, based on 
calculations of unified atmosphere/ wind models with the 
ISA- WIND code (de Koter et al. 1997), confirm this to be 
generally true to good accuracy for the important part of 
the wind (from about the sonic point up to a few stellar 
radii) . Therefore, the complex physics of ion populations 
can be ignored in the investigation of the various effects 
of rotation on the stellar wind. 

However, there is one notable exception which is called 
the bi-stability jump. This refers to a jump in Vqc around 
the temperature of 21000 K where Woo/i'esc climbs from 
1.3 (lower temperature) to 2.6 (higher temperature). This 
jump was observed by Lamers et al. (1995) in B super- 
giants. It is linked to a shift in ionization states of Fe in the 
lower part of the wind. The line driving in the lower part of 
the wind is dominated by iron. Below about T c ff ~ 25000 
K, Fe IV recombines to Fe in and since Fe in is a more 
efficient line driver than Fe iv, the wind structure changes 
dramatically (Vink et al. 1999 ). 
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The bi-stability jump may also occur in the tempera- 
ture difference between the pole and the equator of a fast 
rotating star with gravity darkening. Therefore two sets of 
force multiplier parameters for the wind will be adopted 
to reflect the sudden change in ionization states (a high- 
temperature set for the pole and a low-temperature set for 
the equator). 

3. The physics of rotation 

Intuitively it is clear that rotation has an effect on the 
shape of a star and on the motion of the gas in the wind. 
An additional effect that the rotation of a star can have 
is the darkening of the equatorial regions of a star via the 
von Zcipcl effect. These effects will modify the wind of a 
star considerably and we will incorporate them in the line 
driven wind theory for rotating stars. 

3.1. The shape of a rotating star 

The shape of a uniformly rotating star with all its mass 
concentrated in the core is determined by the equipoten- 
tial surfaces of the potential in a rotating frame (Roche 
model) : 

$(s,M) = - + sin 2 (6) (1) 

x 2 

where x = r/R eq , 9 is the co-latitude (9 = at the pole) 
and to — w ro t,cq/ w crit- is of course independent of the 
longitude (p. Note that l + 0.5w 2 = R eq / 'i? po ic- So the max- 
imum oblateness is R eq /R po i c = 3/2. The critical (break- 
up) velocity is defined in this paper as: 

v 2 crit = GMf/R cq (2) 

where M° ff is the effective mass of the star, 

Mf = M*(l - T e ) (3) 

which is the reduced mass due to radiation pressure by 
electron scattering (see below). 

The surface of a rotating star is implicitly given by 
Eq. [I]. Solving this expression for x(9) is equivalent to the 
solution of a 3rd degree polynomial, which results in 



x{9) = 2 



V2- 



V3 lo sin 9 



■ arcsm 




Because of the difference in radius 

Rr 



R 



= l + 0.5w" 



pole 



(4) 



(5) 



and the difference in rotational velocity between pole and 
equator, the radial escape velocity at the equatorial region 
is smaller than at the pole 



1 - 0.5lo 2 



pole 



i . • > ' sr (6) 

1 + 0.5liH 

We can then roughly estimate the effect of rotation on 
as this value has an almost linear dependence on v csc (see 
later in $t|). 
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Fig. 1. Gravity darkening for different rotation rates lo as 
a function of co-latitude 8 

3.2. Von Zeipel gravity darkening 



The von Zeipel theorem (1924) for distorted stars states 
that the radiative flux from a point on the star is propor- 
tional to the local effective gravity: 



F(9,(b)=a B T^(9) rX g oS (9) 



(7) 



with <tb the Boltzmann constant. As g e s = — V$, we can 
write down the flux of a rotating star as a function of the 
co-latitude 9. This ( lengthy) expression for the flux as a 
function of 9 is plotted in Fig. [j] for a few different values 

of CO. 



3.3. The equation of motion of a line driven wind of a 
rotating star 

The rotation of the star can induce a 9 component of the 
flow of the matter in the wind. In this study a possible 
9 component is neglected as was discussed in §|[ but the 
motion in the longitudinal <f> direction must still be taken 
into account. 

In the absence of forces that can exert a torque on 
the wind, the equation of motion for is given by the 
conservation of angular momentum. Note however that 
there could be a torque from the line forces themselves, 
which is dependent on the velocity gradient of the gas. We 
assume that there are no external torques, so conservation 
of angular momentum gives: 



v<f,(r) = v Iot (Ri,)^y- 



(8) 



The equation of motion for the radial direction be- 
comes: 



dv GM+ 1 dp 



Ri 



dr 



p dr r 6 



(9) 



where v is the radial velocity and M* is the stellar mass. 
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The conservation of mass for a non-spherical sectorial 
wind can be written as 



4irr 2 p(9,r)v(9,r) = F m (0) 



(10) 



where F m (9) is the "local mass loss rate," i.e. the total 
mass loss rate if the solution for this latitude where valid 
for a spherical star. The total mass loss rate from the star 
is 



M 



I 



r/2 



F m (6) sin(9)d9 



(11) 



The equation of motion together with the conservation 
of mass governs the dynamics of the wind given the equa- 
tion of state, p = a 2 p (a is the isothermal sound speed), 
and the radiative acceleration <7 rac i . Using the conservation 
of mass, the pressure term can be rewritten as: 



dv 
dr 



(12) 



1 dp _ 1 da 2 2a? 
p dr p dr r 

As in the CAK-theory the temperature structure is 
specified a priori as a function of r. (The results depend 
only very weakly on the chosen temperature structure; see 
Pauldrachet al. |1986| ). 

The radiative acceleration consists of two components. 
The continuum component due to the electron scattering, 
g e , and a second component, <?l, due to line scattering and 
absorption processes. The continuum acceleration is given 
in terms of T c : 



a n F 



fjc 



= r, 



GAL 



(13) 



where a e is the electron opacity, L is the stellar luminos- 
ity and F is the radiation flux. The contributions by other 
continuum opacities are small. Nevertheless, they were in- 
cluded in the calculation of the force multipliers by Vink 
et al. (1999) used in our bi-stable wind models (Sect. 6). 
For a homogeneous spherical star, T c is given byjj 



4ttGM* 



(14) 



In the more general case of a non-homogeneous, non- 
spherical star the continuum acceleration can be defined 
as a correction to the 'classical' continuum acceleration of 
Eq. M 



T'=D C 



<7 C L* 

j 4ttG?M* 



(15) 



where D c is the continuum correction factor which is given 
by: 



D r 



L * 



I(9,<t>)dn 



(16) 



disk 



1 Note: we define F e as the ratio between the continuum force 
and gravity and not as the ratio between continuum force and 
the critical radiation force for a rotating star 



The line acceleration gi,(r,v,v') has a more complicated 
form, since it is a function of distance r as well as the ve- 
locity v and velocity gradient dv/dr. Combining the equa- 
tion of motion (Eq. ^ with the rewritten pressure and 
continuum acceleration terms and multiplying by r 2 , one 
finds: 



Fg(r,v,v') 



r 2 v 



dv 
dr 



1? 



+GM*(1 - T' c ) - 2a 2 r 



r 2 gh 







(17) 



This equation is valid for each co-latitude 9. Follow- 
ing CAK, the critical point of this equation is found by 
imposing the singularity condition (where the subscript c 
indicates values at the critical point): 



dFe 
dv 







and the regularity condition: 



dFe 
dr 



jt dFe 
dv 



= 



(18) 



(19) 



Note that the critical point is not the sonic point r s 
( v ( r s) — a) due to the fact that there is an additional de- 
pendence on v' in the line acceleration <?l (Abbott 1980; 
ISW Chapt 3.3). 

3.4- The radiative line forces 

The line forces are described within the framework of CAK 
in the Sobolev approximation. This means that the intris- 
tic line absorption profiles are considered to be infinitely 
sharp. Then the line force becomes a function of local 
properties of the wind only. The acceleration due to an 
ensemble of lines is given by the summation of the contri- 
butions of all indivi dual lines with rest frequencies vi. 
is given by (Castor 1974 ): 



9h 



Et 



F 



1 



-pdfl 



(20) 



where ni is the absorption coefficient per gram for the Ith 
line, p — cos 9, I Vl the intensity at the rest frequency vi 
and t Vi is the Sobolev optical depth of Zth line, defined 
as r n — K iP~^(^;)- The integration is performed over the 
complete visible disk of the star. 

Following CAK, the line acceleration can be rewritten 
in terms of the force multiplier as a function of the optical 
depth parameter t. Where t is defined as: 



t 



ref 



dr 
dv 



(21) 



with vth being the thermal velocity of the protons and a™ 1 
is some reference value for the electron scattering <rj cf = 
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0.325 cm 2 g 1 (see ISW, Chapt. 8). The line acceleration 
in terms of the continuum acceleration is give by: 



g L = g c M(t) 



(22) 



If the star is assumed to be a point source, M(t) can 
be approximated by a simple power law parameterization 
(CAK, Abbott 1982): 



(23) 



where the (n e /W) s accounts for the effect of the electron 
density on the ionization balance in the wind. Here n e is 
the electron density in units of 10 11 cm" 3 and W is the 
geometric dilution factor. 

Friend & Abbott (1986) and Pauldrach et al. (1986) 
showed the importance of the finite disk correction on the 
line force in case of an extended source. This finite disk 
correction factor Dfd is given by (CAK): 



4-1. Simplified solutions for non-rotating star 

4.1.1. The point source approximation 

If the star is considered to be a point soure the integral of 
the line acceleration in Eq. collapses to one point. In 
this case Dfd disappears from Eq. ^6] and the line accel- 
eration becomes a simple function g^ = ^(r 2 v^) a . The 
solution of Eq. [D] was found by CAK and can easily be 
explained by neglecting the gas pressure terms ( a 2 ). In 
this case Eq. [IT] reduces to 



r 2 v 



dv 
dr 



GM*(l-r e ) + C 



r 2 v 



(28) 



where C is a constant containing the mass loss rate. The 
equation is solved by imposing uniqueness of the solution 
(Kudritzki et al. 1989) which then fixes the value of the 
constant and thus M. The solution for the mass loss and 
velocity law is given by 



D id = 



M(t) 



AL 



point 



(*) 



l 

N 



L(9, 



0+g) 

1 + CTfi 2 



where 



r dv 
v dr 



1 



and N 



L+R+ 



fuKl (24) 



(25) 



This expression is completely general: the actual shape of 
the star and its intensity distribution enter through the 
integration domain and the 1(9, (f>) dependence. However, 
here we assume the intensity / to be locally isotropic, i.e. 
we neglect limb darkening. We have tested the effect of 
limbdarkening. The inclusion of limbdarkening results in 
a small 9- independent correction on Dfd. 

The line acceleration g^ including all the neccesary 
correction factors is now given by: 



Combining Eqs. [l7| and |26| gives the full equation of motion 
for the sectorial wind of a rotating star 



dv 



dv 



2 / ^ r 



GM i ,(l-T c )-2a 2 r- 



vl 



Rl 



Gtt) 



(vthF m y a D 



id 



(-Y 



r 



4. Solutions of the equation of motion 

For line driven winds of non-rotating stars solutions of the 
equation of motion have been retrieved by e.g. CAK and 
Pauldrach et al. (1986). These will be presented to serve 
as an illustration of the solution of the equation of motion 
and as a basis to interpret the more complicated results 
from the full equation of motion ( p7| ) . 



Al 



47r / cr e \ / 1 — a 



o- e v t h V47T/ V a 

i 

LA 3 



(fca)< 



rn.y 



{GM*(i-r e )}- 



v(r) = Voo [ 1 



with 



a / a 2GM is (l-T) 

T^ Vcsc ~ V (i - a) 



(29) 
(30) 

(31) 



This simplified solution is equal to the full CAK solution, 
in the limit of small sound speed, a « v csc (see also ISW 
Chapt 8). 

4.1.2. Simple finite disk correction 

In case of a homogeneous spherical star the finite disk cor- 
rection factor Dfd (Eq. |J) can be calculated analytically 
(CAK): 



D 



id 



(1-M*) 



(1 



(1 



/idfi 



(l- M 2 )(a + l)a(l + ( 7)« 



(32) 



Including the finite disk correction results in an increase 
in Voo, a decrease in M and a modification of the sim- 
ple scaling laws that were found in the original CAK 



approach (see Friend & Abbott 1986; Pauldrach et al 



1986). The decrease of the mass loss rate compared to 
the point source case is due to the decrease in g^ close to 
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Fig. 2. Geometry for the calculation of the correction fac- 
tors -Dfd and D c . A point on the ray with co-latitude 
sees a different limb angle for every (j) (not drawn; it is 
the angle around the line from stellar center to observer, 
not to be confused with the stellar longitude). The inte- 
grals in Dfd and D c are over < < 27r and < 9 < 0. 
These are rewritten to integrals over the star centered an- 
gles 9' (0 < 9' < 0') and <f> = 4>. The stellar radius and 
flux are given as a function of 9" , the co-latitude, which 
can be calculated for every 9' and </>: 9" = 9" (9', </>). 

the star where M is determined. Close to the stellar pho- 
tosphere the finite disk correction is smaller than one, viz. 
Dfd{r — R*) = 1/(1 + ot) < 1. The accompanying increase 
in of typically a factor of two is due to two effects: (1) 
a reduction of M results in a smaller amount of material 
to be accelerated and (2) far from the photosphere, the 
correction factor Dfd becomes larger than 1. The result- 
ing dependence between and v csc is approximated by 
Friend feAbbott (1986): 



3 / v o 



2.2 



l-aj V^km/s 



0.2 



(33) 



They also found a relation for the modified mass loss rate: 

-0.3 



M ss 0.5 Mcak 



10 3 km/s 



(34) 



4-2. Solution of the equation of motion for the wind of a 
rotating star 

The solution of Eq. |2?j gives the mass loss rate, or rather 
the local mass loss rate F m {9), and velocity structure of 
the wind. The solution is complicated however by the pres- 
ence of D^. This is an integral of a velocity-dependent 
function times the surface intensity over the visible section 
of the star. The appearance of the star varies throughout 
the wind in intensity distribution as well as in shape, as 
the star is no longer spherical and T c ff is a function of lati- 
tude. Since the analytic solution would be cumbersome, a 



numerical approach of the solution is chosen. This is done 
with the analytic solutions of simple, i.e. the non-rotating, 
models in mind for comparison. The numerical solution of 
Eq. ^ is relatively straightforward in case the function 
Dfd is a given function of r. Therefore we solve Eq. [27] 
with Dfd(r) in an iterative way as follows: 



3. Eq 



A (3 law for v(r) is assumed, viz. v(r) = (1 — R+/r) - 5 . 
Dfd(r) is calculated using this velocity law. Note that 
Dfd is not dependent on the actual value of v(r) but 
on the velocity gradient. 

27 is solved using this correction factor Dfd(r), 



obtaining a new velocity law v(r). 
4. Step 2-3 are repeated until convergence is reached. 

Typically three iterations are sufficient, since further 
iterations changed the obtained values by less than 0.5 
%. We checked that the solutions converged to the same 
value, independent of the starting value. 

4-3. The calculation of Dfd (r) and the continuum 
correction factor D c 

The calculation of the correction factor Dfd from the ve- 
locity law (step 1. in the scheme described above) is per- 
formed by numerical evaluation of Eq. E4l This is a non- 
trivial task, since both the integrand (Eq. |7j) and the shape 
of the visible "disk" (Eq. ||) are dependent on r. 

The main parameters needed to obtain Dfd for a fixed 
co-latitude 9 are: the stellar surface R(9) and the sur- 
face temperature T(9), both given as a function of the 
co-latitude 9, as well as the velocity law v(r). It is conve- 
nient to change the integration variables in Eq. ^Jfrom the 
wind centered coordinates 9 and <p to the star centered co- 
ordinates 9' and 4>. For the definition of these coordinates, 
see Fig. |2|. Note that here 4> is not the stellar longitude 
anymore, but the angle of rotation around the line from a 
point (r, 9) in the wind to the stellar center. 

1 f*" f e( -^ ( (1 + a) \ a 
N J J \1 + cr cos 2 9 J 

(35) 



1 



2tt r e'(<p) 



, (1 + " } 2 Vcosflsinfl^dfl'c 
1 + a cos 2 9 J dO' 



N Jo jo 

where I is the frequency-integrated intensity 
I = 1(9', 0) oc T{9" f (36) 
The angles are related to one another via 
r-R(6',(j>) cos 9' 



cosy = 



(37) 



yj(r - R{9' : (j)) cos 9') 2 + (R(9' , (j)) sin 9') 2 
with the stellar surface given by R(9' , (f>) = R(9"), and 
cos 9" — sin 9'coscf> sin oj + cos 9' cos w (38) 
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Equation |3^ can be integrated, considering that the limb 
angle &'{<p) can be determined numerically, yielding the 
finite disk correction for the line acceleration of an arbi- 
trary shaped star. 

The correction factor D c for the electron acceleration 
(defined in Eq. Hq) similarly becomes: 



1 



2tt ,-e' 



d6 



D c = — I I I cos9 sin 9 — d9'c 
N J Jo d0> 



(39) 



This correction is also included in the solution of the equa- 
tion of motion. 



4-4- Solving the equation of motion 

To solve the equation of motion it is neccesary to deter- 
mine the conditions at the critical point: r c , v c and v' c , or, 
equivalently, r c , v c and vo — v(R+), as well as the local 
mass loss rate F m . Apart from Eqs. |l^ and [l^ and the 
equation of motion (|27]) itself the additional constraint 
for the continuum optical depth r G at the stellar radius is 
neccesary (Pauldrach et al. 1986) to uniquely determine 
these quantities: 



pcr dr 



(40) 



If the condition for vo is used one may adopt any small 
value, say Vo = 0.1 km/s, since even a large error in vo will 
only yield a relatively small error in Voq and F m . This is 
due to the fact that at the surface of the star v(r) and thus 
also p vary exponentially with scale height H = 



a 2 Rj 



« 

. This means that the radius where Eq. |40| is fulfilled 
is very close to R*. This has been checked a posteriori. In 
our calculations we have adopted this boundary condition 
for vo in the solution of the full momentum equation. 

A first guess for the local mass loss rate obtained from 
the CAK solution is used to integrate the equation of mo- 
tion from the stellar surface outward. This gives an ap- 
proximate value for the critical point where r 2 v^ = 
if F m (guess) < F m (solution). This approximate value for 
the critical point can then be used to consecutively im- 
prove the guessed value of F m using Eqs. |l8| and This 
iterative procedure is terminated when the local mass loss 
rate converges. The calculation is then completed by in- 
tegrating past the critical point on the solution which ex- 
tends to infinity (see CAK or Abbott 1980 for a description 
of the topology of the solutions). To deal with a numerical 
singularity at the sonic point the equations are solved for 
r(v) rather than for v(r). 

5. Application to B[e] winds 

5.1. A typical B[e] supergiant 

We have calculated models for rotating B[e] stars using the 
above described method. For the first models we adopted 






0.3 (pole) 
0.3 (equator) 



r/R 




Fig. 3. The velocity law v(r) at the pole and equator of 
B[e] wind models for different rotation rates. Top panel: 
oj = 0.3, bottom panel: uj = 0.6. The crosses indicate the 
location of the critical point of the wind. Note that the 
equatorial curves in the rotating models start at a radius 
r/Rpoic > 1 due to the oblateness of the star (see text for 
stellar parameters). 



the following stellar parameters, which are typical for B[c] 
stars: L* = 1O 55 L , M* = 4OM , R^ pole = 47i? Q , 
T e g ~ 20 000 and solar abundances. The effective temper- 
ature is defined as T e g = (L^/o-bS) - 25 where S is the total 
surface of the distorted star. The following force multiplier 
parameters from Pauldrach et al. ( |1986 ) were adopted 
to describe the line acceleration: a = 0.565, k — 0.32, 
5 = 0.02. An overview of the results for this generic B[e] 
star rotating at ui = 0, 0.3, 0.4, 0.5 and 0.6 times the crit- 
ical rotation speed is given in Table 1. The calculation of 
models with higher values of u> runs into numerical prob- 
lems, because in the iteration, the approximate critical 
point can be below the sonic point. The table shows for 
different co-latitudes 9 and different rotation velocities u> 
the values of Voo, the local mass loss rate F m {9) and the 
value of the velocity law parameter /3, obtained by fitting 
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Table 1. Properties of a typical rotating B[e] star model 1 
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1.40 


2.35 


0.63 


























.3 


1.43 


2.67 


0.63 


1.40 


2.64 


0.63 


1.35 


2.53 


0.63 


1.30 


2.42 


0.63 


1.29 


2.36 


0.64 


.4 


1.44 


2.95 


0.63 


1.40 


2.88 


0.62 


1.31 


2.70 


0.62 


1.23 


2.46 


0.63 


1.20 


2.34 


0.64 


.5 


1.46 


3.33 


0.63 


1.40 


3.23 


0.62 


1.27 


2.95 


0.62 


1.14 


2.53 


0.64 


1.11 


2.27 


0.66 


.6 


1.48 


3.83 


0.63 


1.40 


3.70 


0.62 


1.21 


3.29 


0.60 


1.04 


2.61 


0.64 


1.01 


2.11 


0.70 



(1) The adopted stellar parameters are: 

T cft = 20000 K, U = 1O 5 ' 5 L , Mi, = 4OM , = 47i? , a = 0.565, k = 0.32, 8 = 0.02, solar abundances. 

(2) : the terminal velocity is in 10 3 km/s 

(3) : the mass flux F m in 10" 6 M y _1 

(4) : the velocity law parameter j3 is obtained from a nonlinear fit of v(r) from 1.1 J?* to 10 R*. 



Table 2. Properties of rotating B^] 1 models: with different values of L+ 





pole ppole eg jpeq Peg Jl> 

If) 3 lem e^ 1 ID -8 M^v~ 1 ID 3 km 1 D~ 8 A/f^ir - 1 ID -8 M^\r~ 1 
_LU Kill b lu ivi y -LU JVIIl fa 1U IvlQy 1U IVlQy 


1Q 4.b 


2.84 3.01 1 3.01 
0.3 2.91 3.37 2.60 3.04 1.01 3.18 
0.4 2.97 3.67 2.42 3.03 1.01 3.32 
0.5 3.04 4.08 2.22 2.98 1.00 3.50 
U.O o.ll 4.0Z z.uu z.oz u.yo o.ll 


u 

L 


poie ppole e 9 peg _Peq_ yw- 

10 3 km s" 1 10 -7 A/ y _1 10 3 km s" 1 10" 7 Moy -1 10" 7 Moy" 1 


10 b.u 


2.06 2.53 1 2.53 
0.3 2.11 2.84 1.88 2.55 1.00 2.67 
0.4 2.15 3.11 1.75 2.54 1.00 2.79 
0.5 2.19 3.46 1.61 2.49 0.98 2.94 
0.6 2.24 3.93 1.46 2.35 0.91 3.15 


U 
L 


w v pol f F£ o!e v eq f F% -^a- M 

mf m inj m p po le 

10 3 km s" 1 10" 6 Mey- 1 10 3 km s" 1 10" 6 Moy" 1 10" 6 Moy" 1 


10 s0 


1.40 2.35 1 2.35 
0.3 1.43 2.67 1.29 2.36 0.98 2.47 
0.4 1.44 2.95 1.20 2.34 0.95 2.56 
0.5 1.46 3.33 1.11 2.27 0.90 2.71 
0.6 1.48 3.83 1.01 2.11 0.80 2.89 


u 


u v p ° l ; fT !e v eq f F% M 

10 3 km s" 1 10~ 5 Mey" 1 10 3 km s _1 10" 5 Moy -1 10" 5 M y _1 


10 b.u 


0.65 3.93 1 3.93 
0.3 0.62 4.97 0.60 3.80 0.78 4.25 
0.4 0.59 6.00 0.58 3.61 0.62 4.54 
0.5 0.55 7.79 0.56 3.24 0.41 5.01 
0.6 0.49 11.1 0.56 2.64 0.21 5.76 



(1): The adopted stellar parameters are: M t = 4OM ; R+ = 47R Q ; T eS = 20000 K. 



the calculated velocity structure from 1.1 R* to 10 i?* to 
a /3-law. We see that the value of j3 does not change much. 

Fig. U shows the velocity structure of the wind and the 
location of the critical point for u> = 0, 0.3 and 0.6. Note 
that the critical point is very close the star. This is true 
for all finite disk corrected solutions. Figure [| shows the 
behaviour of versus F m for various rotational speeds 
and latitudes. As expected, at the equator is smaller 



than at the pole because Voo scales roughly with v csc and 
due to the smaller value of u csc at the equator. 

Figure || shows the latitude dependence of the local 
mass-loss rate and of the terminal velocity from equa- 
tor to pole, for different rotational speeds, indicated by 
to. At the pole (9 = 0) Uqo is about constant for various 
rotational speeds. The mass- loss rate at the equator de- 
creases slightly with increasing rotation. The increase of 
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Table 3. Properties of rotating B[e] 1 models: with different values of M* 
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1.76 


1.66 


1.56 


0.99 




1.55 




0.4 


1.89 




1.92 


1.56 


1.55 


0.98 




1.59 




0.5 


1.91 
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1.51 
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1.95 
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(1): The adopted stellar parameters are: L* = 10 5 5 L Q ; R* = 47i? Q ; T off = 20000 K. 
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Fig. 4. The latitude dependence of the mass-loss rate and 
of a typical B[e] star, for various rotational speeds. 
(See text for stellar parameters). 
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the mass los s rate at the equator as calcu lated by Paul- 
drach et al. ( |1986| ) and Friend & Abbott (|l986| ) is offset 
by a decrease in mass loss rate due to the smaller radiative 
flux at the equator due to the von Zeipel effect. The mass 
loss rate at the pole increases strongly with increasing ro- 
tation rate due to the increase in luminosity at the pole. 
This is because the total luminosity of the star is fixed 
and the smaller luminosity at the equator must therefore 
be compensated by a higher luminosity at the pole. 

5.2. The overall density properties 

The overall effect of rotation is to increase the mass loss 
rate of a star. This can be seen in Fig. [3] where the overall 
mass-loss rate (i.e. F m integrated over the stellar surface) 
is plotted vs. the rotational velocity. The mass-loss rate 
of the models with 4.5 < logL*/L0 < 5.5 vary with w as 
1 + 0.64W 2 . 

Figure ^ shows the resulting density contrast (far 
from the star) between pole and equator ( p eq / p po \c = 
Fm v Jx lc /Fm lCv( ^)- The densities at pole and equator are 
essentially the same: the smaller Vqo at the equator is offset 
by the larger mass loss rate F m at the pole. 

5.3. Varying L+ 

We have also calculated models for different luminosities 
(see Table 2), and different stellar masses (see Table 3). 
The mass-loss rate for the different luminosities is also 
plotted in Fig. |^. The density contrast for various lumi- 
nosities is plotted in Fig. [| We see that the density at 
the pole compared to the equator increases for increas- 
ing luminosity and the effect on the mass-loss rate is only 
visible for the highest luminosity. This is an effect of the 
continuum acceleration D c . Fig. 7 illustrates the effect of 
the continuum acceleration on the forces in the wind for 
a star with high luminosity (r e = 0.76). The net effect 
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Fig. 6. The density contrast p c(i /p polc 
of log L^/L Q . 



for various values 



/ r f r 
/ / / 




Fig. 7. Vector plot of the gravity (F g ) and the continuum 
radiation force (F c ). The vectors indicate strength and 
direction of r 2 (F g + F c ) on an arbitrary scale for a star 
with r e = 0.76 and uj = 0.6. Note that line forces and 
centrifugal forces are not included. 



of the continuum radiation pressure and the gravity, i.e. 
GM*(1 — T' e )/r 2 , is an inward force in the equatorial re- 
gion and an outward force in the polar region where the 
flux is higher. 

For small luminosity the variation in T' c through D c is 
not important because T c itself is small. For high luminos- 
ity L*/(GM+) increases and thus the effect of D c on T' c 
becomes noticeable (this is with all other properties of the 
star fixed) . In Eq. we see that a larger value of D c at 
the pole causes a relatively larger decrease (through T' e ) 
in v esc and thus in Voo, whereas the mass loss rate at the 
equator will increase less because of the smaller D c . This 
is only a qualitative explanation and the actual depen- 



dence of Peq/Ppoie does not follow the relation suggested 
by Eqs. 29 and 31. We have found that the actual values 
of p eq / Ppoio could be approximated quite well (within 10 
percent) with the following relation: 



Peq/ Ppolc 



l-r' c {9 = 0,r/R cq = 1) 



l-T' e (6 = ir/2,r/R. 



eq 







1.5 



(41) 



where the r m i n « 1.2 is the radius where D c reaches its 
minimum value. 

So we see that although rotation alone modifies the 
wind structure considerably, it cannot be responsible for 
the density contrast between the equatorial and polar 
wind observed in B[e] stars, in case only radial effects are 
considered. We conclude that at least one other physical 
effect must be responsible for the formation of disks of ro- 
tating B[e] stars. The two most likely additional effects are 
the bi-stability and the flow of the wind material towards 
the equator. 

6. Rotationally induced bi-stability models 

In the section above we have shown that rotation alone 
cannot explain the observed large density contrast be- 
tween the poles and the equator of B[e] stars. One of the 
possible additional mechanisms to enhance this contrast 
is the bi-stability jump. 

To investigate whether the bi-stability jump observed 
in normal B-supergiants can explain the formation of disks 
around B[e] stars, we have calculated bi-stable wind mod- 
els for a fast rotating typical B[e] supergiant. Spherical, 
NLTE wind models for normal B supergiants have been 
calculated by Vink et al. (1999). From these models CAK- 
like force multiplier parameters have been derived using a 
Monte-Carlo method to simulate photon-gas interactions. 
These models show the occurence of a bi-stability jump 
around T cS ~ 25000 K. 

For two models on either side of the bi-stability jump 
(T cff = 17 500 and 30 000 K), the CAK parameters were 
determined by fitting a power law from about the sonic 
point to about 0.5voo (from t = 10~ 2 to t = 10~ 4 in 
optical depth, t is defined in Eq. ^l|) to the calculated 
force multiplier values. The resulting values for fc, a and 
5 are listed in Table |^. Note that the <5 parameter is taken 
equal to zero as it was not possible to extract it explicitly 
from the models. Its effect is hidden in the constant fc. As 
expected, the values for the force multiplier parameters 
are quite different for the two cases, because the ionization 
has changed dramatically over the bi-stability jump. 

The predicted bi-stability jump occurs around T c s ~ 
25 kK. Therefore we calculated a model for a rotating 
B[e] star with the following properties: T e g = 25000 K, 
L± = 1O 5 L , M± = 17.5-Mq and solar abundances. This 
implies that the pole of the rapidly rotating star will be 
hotter than 25 kK and the equator will be cooler than 25 
kK. 
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Table 4. 

wind 

T c ft 



The force multiplier parameters for a bi-stable 
k a 8 



17500 K 
30000 K 



0.57 0.45 0.0 
0.06 0.65 0.0 
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Fig. 9. The mass loss rate F m of a star with a bi-stability 
jump as a function of co-latitude. The dashed part in- 
dicates the location of the bi-stability jump. Notice the 
steep jump in mass loss rate near 6 — 60°. (See text for 
the adopted stellar parameters) 



Fig. 8. The terminal velocity Voo for a star with a bi- 
stability jump as a function of co-latitude. The dashed 
part indicates the location of the bi-stability jump. Notice 
the drastic drop by about a factor 3 around 6>=60°, due 
to the bi-stability jump. (See text for the adopted stellar 
parameters) 



For the pole, the force multipliers k and a from the 
hot Monte Carlo model of 30 000 K were used, whereas 
for the equator, k and a from the cool model of 17 500 K 
were used. 

The resulting Voo(9) and mass loss rates F m (9) are 
plotted in Figs. || and || for oj = 0.6. Clearly visible is 
the drastic decrease in and the drastic increase of F m 
towards the equator. This occurs around a co-latitude 9 
(angle from the pole) of about 60 degrees, where T e s ~ 
25 000 K. (The precise location of the jump depends on 
the effective temperature of the star) . Since both F m and 
Voo have a equator/pole contrast of about a factor of three, 
the resulting density contrast is about a factor of 10 in a 
disklike region with a half opening angle of 30 degrees. 

Our calculations show that the density contrast in the 
wind between pole and equator due to the bi-stability 
jump is significantly larger than without the bi-stability 
jump. Yet, the calculated value of about a factor 10 is 
not sufficient to explain the observed density contrast of 
a factor ~ 100. So another mechanism is needed to en- 
hance the density contrast even further. This is most likely 
the wind compression mechanism (Bjorkman & Cassinelli 
1993; Bjorkman 1998). 



7. Summary and discussion 

We have modified the line driven wind theory for rotating 
stars by including the effects of oblateness and gravity 
darkening as well as rotational terms in the equation of 
motion. We considered a sectorial wind, i.e. we neglected 
the effect of the motion of the wind towards the equator or 
towards the pole. This assumption is justified close to the 
stellar surface below the critical point of the momentum 
equation. This implies that our method will predict about 
the correct distibution of mass-loss rates from the star as 
a function of stellar latitude, but it may not be accurate 
enough to predict the velocity and density distribution 
further away from the star, if motions in the 9 direction 
become important. 

The equation of motion was solved using an iterative 
numerical scheme, that includes the conservation of an- 
gular momentum and the correction factors to the radia- 
tive acceleration by lines and by the continuum, due to 
the non-spherical shape of the star and due to the lati- 
tude dependence of the radiative flux. This method was 
applied to study its possible effects on the formation of 
disks around fast rotating B[e] supergiants. 

The models with constant force multiplier parameters 
k, a and S show a decrease of both the mass loss rate and 
the terminal velocity from pole to equator. This is mainly 
due to two effects: the reduction of the escape velocity 
from pole to equator, resulting in a higher terminal veloc- 
ity at the pole, and the reduction of the radiative flux from 
the pole to the equator due to gravity darkening, which 
results in a decrease in mass loss rate at the equator. 

For a star with a fixed luminosity and a fixed polar 
radius, the temperature at the pole increases and the tem- 
perature at the equator decreases with increasing rotation 
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rate. The terminal velocity of the wind from the poles is al- 
most independent of i> rot but at the equator Voo decreases 
with increasing rotation rate. The mass-loss rate at the 
pole increases with increasing v rot , due to the increase in 
T e ff at the poles, but the mass loss rate from the equator 
is almost independent of the rotation. The combination of 
these effects alone produce a density contrast between the 
polar and the equatorial wind of a factor p e q/p po ie — 1, ex- 
cept for fast rotating stars with high luminosity (in which 
case Poq/Ppoic < 1)- 

Our results agree quantitatively with those obtained 
by Maeder (|1999|) for the latitudinal dependence of F m . 
Our predictions of the global mass loss rates are also in 
general agreement with those of Maeder for models with 
luminosities far from the Eddington limit. For high lumi- 
nosities we find a strong polar outflow, whereas Maeder's 
models predict an enhanced equatorial outflow on the ba- 
sis of the original CAK scaling laws. This difference is due 
to the fact that we have included an improved description 
of the radiation pressure in a rotating star. 

The difference between the polar wind and the equa- 
torial wind is strongly enhanced when the bi-stability of 
radiation driven winds between T c ff ~ 20 000 and 30 000 
K is taken into account. In this case the force multipliers 
k and a of a rapidly rotating star change drastically with 
stellar latitude if the pole is hotter than 25 000 K and 
the equator is cooler. In § [| we have applied the newly 
calculated force multipliers from Vink et al. (1999) above 
and below the bi-stability jump to the models of rotating 
B[e] supergiants. One might argue that it is not allowed 
to apply the force multipliers of (bi-stable) spherical wind 
models to aspherical winds of rotating stars. However, we 
have shown that a disk formed through bi-stability will be 
sufficiently thick (typically about 30 degrees for a star with 
an effective temperature equal to about 25 000) to make 
these spherical models reasonable approximations for the 
conditions in the wind close to the star, where the mass 
loss rate is determined. We find that rotationally induced 
bi-stability models of B[e] stars reach a density contrast of 
about a factor 10 between the dense equatorial wind and 
the less dense polar wind. This is le ss tha n the factor 10 2 
that is observed (see e.g. Bjorkman 1998 ). 

The extra density increase is most likely due to the 
wind compression. The gas that leaves the photosphere 
from a rotating wind, will follow an orbit in a tilted plane 
defined by the local rotation vector and the center of the 
star. If the rotational velocity is large or the wind velocity 
is small, this orbit will cross the equatorial plane where 
the streamlines of the wind from different stellar latitudes 
cross. The resulting shock will compress the gas into a 
thin outflowing disk with an opening angle of only a few 
degrees (Bjorkman & Cassinelli 1993; ISW Chapt 11). We 
have not considered this wind compression in our model. 
However, it is likely to occur in rotationally induced bi- 
stable winds, because the wind velocity in the equatorial 



plane is about a factor three smaller than from the poles, 
thus facilitating the wind compression. 

Owocki & Cranmer (1994) have argued that the flow 
of the wind towards the equatorial plane, predicted in 
the WCD-theory, may be offset by a 0-component of the 
line acceleration towards the polar regions. This "disk in- 
hibition" mechanism operates in wind models with con- 
stant force multipliers k and a. However, it is not clear 
that this mechanism is sufficiently strong to overcome the 
wind compression in the RIB-model. This is because k and 
a change with latitude in the RIB-model and there is a 
strong density gradient in the wind from the equator to 
the poles (Owocki et al. 1998, Puis et al. 1999). 

The combination of the RIB and the WCD mecha- 
nism may offer the best possibility for explaining the disks 
of B[e] -supergiants. Whether the compression is strong 
enough to explain the observed high density contrast be- 
tween the polar and the equatorial wind, remains to be 
seen. It can be calculated by combining the solutions of 
the wind momentum equation of rotating oblate stars with 
gravity darkening (derived in this paper), and the calcula- 
tion of the resulting trajectories of the wind. The combi- 
nation of the rotation induced bi-stability model and the 
wind compressed disk model is promising for explaining 
the disks of B[e] stars because the RIB-mechanism ex- 
plains the increased mass loss and the small velocity from 
the equatorial regions and the WCD-mechanism explains 
the strong compression of the disk. 
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